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Abstract 

We present an analysis of the gonihedric loop model, a reformulation of the two dimensional 
gonihedric spin model, using two different techniques. First, the usual regular lattice statistical physics 
problem is mapped onto a height model and studied analytically. Second, the gravitational version of 
this loop model is studied via matrix models techniques. Both methods lead to the conclusion that 
the model has Cmatter ~ for all values of the parameters of the model. In this way it is possible to 
understand the absence of a continuous transition. 

1 Introduction 

The gonihedric string model (a type of string model with vanishing tension at tree level) was introduced 
and elaborated further in ,1 . The proponents argued its potential relevance to describe the QCD string 
and is, in any case, an interesting model of random surfaces. The so-called gonihedric spin model provides 
a dual version of the gonihedric string embedded in a (hyper)cubic discrete spacej^lEl- This spin model 
describes the gonihedric string as the interface between positive and negative spins regions. 

The thermodynamical properties of this spin model have been already considered before. The existence 
of a second order phase transition separating the 'ordered' and 'disordered'^ phases at T = Tc would make 
it possible to describe the continuum properties of the discrete gonihedric string. Several papers studied 
the phase diagram of this model .5 , & , but the basic version of the gonihedric spin model exhibited a first 
order phase transition rather than a continuous one. This led the inclusion of a new parameter k 
that regulates the self-avoidance of the 'string'. For k — the original non-self-avoiding string theory is 
recovered, while for increasing k the string becomes more and more self-avoiding. Of course a non-zero 
value for k is probably of physical interest anyway. With the addition of this new parameter the phase 
diagram changed considerably exhibiting second order phase transition for k greater than some Kc- 

At some point the spin model gained interest by itself. Extremely slow approximation to the equi- 
librium in a region below strongly suggested the presence of glassiness in the dynamics of this three 
dimensional model. The model was studied in 0. It revealed indeed the presence of glassiness, when 
K = 0, in its dynamics for T < Tg, where Tg {Tg < Tc) is the glassy temperature. Tg was determined 
with good precision in the model. For k 7^ it was seen that the first order transition present at k = 
persists up to some Kc- From this point on the transition is of second order. Also Tg tends to Tc for k 7^ 
up to a point where it becomes impossible to tell them apart; for big enough k the glassy phase extends 
all the way up to Tc. The dynamics of a family of very similar models (but not including the gonihedric 
model for any value of k) had already been studied in 8. . 

The very geometrical origin of this spin model fine tunes the coupling constants (only nearest neigh- 
bors, next to nearest neighbors, and plaquette-like operators appear) to very concrete values. With the 
addition of the self-avoidance k parameter, a one parameter family of spin models survives. In particular 
the energy of a spin configuration can be seen to be 

E ^ n2 + Kn4 (1) 



^E-mail: espriu@ecm.ub.es 
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•^The quotation marks are appropriate here as this model has a exponentially large degeneration and so what order and 
disorder means is not completely apparent. See 01 
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where 712 is the length of corners of the inter-phase separating plus and minus spins, and is the length 
of the self crossings of this inter-phase, all in the lattice units. This is in fact the energy of the gonihedric 
string, where the flat regions of the string does not contribute to it at all. 

These fine tuned coupling constants are the reason for some crucial features. The Hamiltonian ex- 
hibits a large amount of symmetry; one can for instance flip at once any whole plane of spins (in three 
dimensions) in any direction without changing the energy of the configuration (some of this symmetry 
is lost when k 7^ 0). This feature is thought to be relevant in the appearance of the glassiness 'Tj. The 
appearance of energy barriers along the evolution of the spin system is another determinant characteristic 
of this system which is certainly responsible for very slow dynamics; some metastable states are formed 
where the system gets trapped for a long time. 

In [5j the model was reduced to two dimensions and the dynamics in the two dimensional model 
studied. Eq. remains valid and now 712 corresponds to the number of corners and the number 
of crossings of the boundary line separating plus and minus domains. This two dimensional model can 
actually be represented by a model of loops. We shall term this the gonihedric loop model. 

Although the two dimensional model exhibits large degeneracy and trapped, long lived, metastable 
states, the dependence of the barriers on the linear scale of the system is different and no glassiness 
was found in this case. Only a very slow, power-like, but certainly faster than logarithmic, evolution to 
equilibrium was found. An analytical estimation of the power-like behavior was computed and found to 
be in reasonable accordance with numerical simulations. 

As for the thermodynamical properties, the two dimensional model seems to have no criticality at all. 
It is in fact solvable and trivial for k = (it actually corresponds to one of Baxter's six vertex models pi!|'l 
while there is no analytical solution for other n values. 

In this work we plan to study with different analytical techniques the gonihedric loop model, both on 
its own and coupled to gravity. In this latter version the underlying lattice is dynamical. 

In section|21we briefly recall the basic definitions of the gonihedric loop model. In sectional we study 
the gonihedric loop model in two dimensions in the framework set by Kondev et al. in jjl and by Di 
Francesco and Guitter in a series of papers (see |12| for a review and references therein) . These techniques 
are developed to be used in loop models with no temperature or energy parameter; they correspond to a 
combinatorial problem. With some heuristic arguments we will connect their techniques and results with 
our model. 

In section 0] we couple the gonihedric loop model to gravity. This can be performed using Random 
Matrix Model techniques, more precisely we will write it in terms of a Hermitian two-matrix model. This 
new model is already interesting by itself as it has not been solved exactly before (although some very 
similar models have been, see The peculiar interactions between the two matrices with independent 

coupling constants makes it very hard to solve. Anyway this will provide us an analytical handle on the 
gonihedric loop model on randomly generated lattices. In the limit n 00 one of the matrices in the 
matrix model reduces to a Gaussian variable that can be exactly integrated out. Unfortunately this does 
not make the model directly solvable. To test then the heuristic arguments we have alluded to before we 
study the renormalization group flow of the coupling constants using the methods in 14 . As a byproduct 
we flnd again confirmation of the non-existence of a continuous phase transition in the gonihedric spin 
model through the KPZ relations )15| . 

2 The gonihedric loop model. 

First of all we shall present how to construct the loop configurations from the gonihedric spin configura- 
tions and how to identify the points where energy is located. 

The gonihedric spin model Hamiltonian has the following form in two dimensions. 



When computing the energy we realize that there are only three possible levels of energy for a plaquette. 
In figure |21 we can see all possible configurations up to rotation and Z2 symmetries. As we can see in 
the figure we can associate to each spin configuration the inter-phase configuration just drawing a line 
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(dashed line) that separates plus and minus spins. This inter-phase will be the gonihedric loop and the 
energy will be accumulated in configurations where the loop bends or crosses with itself. 
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Figure 1: All kind of plaquette terms up to symmetries and its corresponding loop segment. Figure 1(a) 



and 1(b) contributes with some cq energy while figure 1(c) contributes with 4k and figure 1(d) with 1 to 
the energy. 



Using the Hamiltonian in Eq. ^ we compute the energy costs to generate bendings and crossings. 
While figure 1(a) and 1(b) cost some bulk energy eo, figure 1(c) costs cq + 4k and 1(d) costs eg + 1 in 
energy units. That means that all the energy will accumulate in regions where the loop is bended or it 
crosses with itself. Then, apart for a configuration independent term proportional to the volume, we get 
an energy per configuration given by equation l|T|l. 

The Hamiltonian ^ thus describes a theory of loops defined on a two-dimensional lattice, modulo 
the usual Z2 degeneracy of Ising-like models. The weight of each loop configuration is given by formula 
This is the Gonihedric loop model. 

Particularly for k ^ the model looks like an ordinary Ising model with some additional interactions. 
This similarity is in fact way too naive. In fact the critical behavior as studied in the references already 
mentioned is quite different. In this work we shall understand why. 

In fact, when placed on a square regular lattice, the above statistical model is trivial for k = (only 
plaquette interactions survive in this case, and we obtain a model quite different from Ising). For k 7^ 
there is no known exact solution of the model. When we place the above Hamiltonian on a randomly 
generated lattice, it will correspond to a model of two dimensional gravity possibly with some matter 
contents. The former case is studied in the next section while the latter is studied in sectional 



3 Relation to height models. 

In a series of works (see |12| for a review) several aspects of loop models in connection with folding and 
meander problems have been studied. There is a mapping from loop models onto the so called height 
models. These height models are going to be interesting to us because they provide a systematic way to 
conjecture the actual value of the central charge c. These results are actually based on a previous work 
by Kondev et. aL jll| where the continuous effective Coulomb gas formulation is studied. 

Let us describe the height models for different classes of loops and recall some of the results and 
formalism by Kondev et. ai. We assign to each loop a weight n — 2cos(7re). The 6 = case {n — 2) 
corresponds to the oriented loop model*. The case e = i, corresponding to a weight n = 1, corresponds 
to non-oriented loops. Let us work first with this e = case, i.e. oriented loops. 

*A weight n = 2 may be understood as two possible equally probable states for each loop, this could be encoded in their 
orientation. 
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3.1 Fully and dense packed loop models on a trivalent regular lattice 



Let us begin by considering the so-called fully packed loop model on a regular trivalent lattice. This 
section is based on J^. Fully packed means that each site is visited once and only once by a loop. We 
will call cfpl the central charge of this model. An example of such a loop configuration is drawn in 
figure m The bicolorability of nodes of the lattice enables us to attach labels on the loop unit segments. 




Figure 2: An example of an oriented loop model configuration on a bicolorable regular trivalent lattice. 
The left one is for the fully-packed model while the right one is for the dense-packed model. The 
corresponding central charges are cfpl and cdpl 

This is equivalent to saying that every loop has an even number of unit segments — then they can be 
labeled using two labels alternatively. We call B the segment joining a black node to a white node 
following the loop orientation, and C the segment joining a white node to a black one following the loop 
orientation. We will give a label also to the unvisited edges, let us say A. In figureElthis label association 
is represented. Once we have defined the different type of unit segments in our theory we will define the 
associated height theory. 

Define a height X variable in the center of every face. This height variable X will differ from one face 
to another, it actually increases by a defined quantity each time we cross a link. This quantity is defined 
following the rules shown in figure |21 Now we just need to ensure that the height variable is globally 
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Figure 3: Rules defining the increments of the height variable when crossing an edge (left) and condition 
for the height increment unit vectors A, B and C when going around a vertex (right). The bottom one 
is present only in the dense-packed loop model. 

defined, that means that the total gained height in any closed path is zero. This amounts in this case 
to the condition A + B + C ~ 0, i.e. two real degrees of freedom remain and so we can describe the 
model with a two-dimensional local field. This implies in the continuum limit a model with two scalar 
fields (the two components of the height field) and so the central charge for the fully packed loop model 
Cfpl = 2. 

A dense packed loop model is one where the condition of each site being visited by a loop is relaxed. 
In the case of dense packed loop model we need to consider a new kind of vertex. This new vertex is 
just the non-visited vertex, the vacancy or impurity vertex. The condition we find after walking around 
a non-visited vertex is A = 0, thus B = — C, and so the associated height model reduces to a one- 
dimensional height variable. In the continuum limit the associated continuum height model contains only 
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one scalar height field and so in the dense case the central charge for the dense-packed loop model is 
CD PL = cfpl — 1 = 1- Notice that the B — —C condition erases all information about the bicolorability 
of the lattice. We can see in figureElthat when B = —C it does not matter what color the nodes are. 

This was treated in general in the works of Kondev, de Gier and Nienhius and used by Di 
Francesco and Guitter 12 . In the general picture a weight n = 2 cos(7re) per loop (e S [0, 1]) is introduced 
and the the model is studied following an effective field theory description. This analysis gives 

CDPL = CFPL -1 = 1- 6- . (3) 

1 — e 

For the non-oriented case, equivalent as we have already said to e = ^ in the general description, we 
can still analyze the height components in the same form. We have no orientation so we cannot use the 
bicolorability information to define in a meaningful way the labels on the segments. Then we get only 
two kind of height increments, say A and B. In figure 01 we can see that there is no way to distinguish 




Figure 4: An example of a non-oriented loop model configuration on a bicolorable regular trivalent lattice. 
Fully-packed loop model on the left and dense-packed loop model on the right. 



different kind of segments on the loop. The bicolorability property combines with the orientable nature 
of the loops to define the B and C increment vectors. The loss of any of the two ingredient forces a 
constraint between B and C. The rules for the height increments are shown in figure |S1 

In the fully packed non-oriented loop model the only vertex we have is again a visited vertex with 
two loop segments and an empty edge. The condition ensuring a meaningful height variable is thus 
A + 2B = 0. Following again the same steps as before we end up with a one dimensional height variable, 
meaning that cfpl — 1 in the non-oriented loop model on a trivalent regular lattice. 

The dense packed non-oriented loop model is again straightforward. The extra condition for the new 
unvisited vertex is 3A = 0. This constraints A and B to zero so the height variable remains constant in the 
whole lattice. This means no scalar field in the continuum theory, and consequently cjjpl = cfpl~^ = 0. 
Note that this coincides with the values we would obtain from equation |(2Jl when applying e = i. 
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Figure 5: Left: In this case the height increment rules involve only two increment vectors, A and B. 
Right: The conditions on this vectors 

These two examples help us to understand how the height models are constructed, and how to extract 
the information we are looking for. 
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3.2 Non-oriented dense loops on tetravalent regular lattice. 



Taking the gonihedric spin model as our starting point we can arrive at a loop model similar to the ones 
we have just described. One difference is that our spin model is placed in a square lattice and so are the 
loops that live on the dual lattice. As already explained above these loops will add energy to the system 
due to bending {AE = 1) and crossing {AE = 4k). 

Let us consider first the case where k = oo (infinite self-avoidance). In this case there is no crossing 
of the loops, so if we forget for a moment about the temperature we have a model of dense loops with 
no self crossings. They are of course not oriented loops since there is no way the spin configuration can 
induce an independent orientation to each loop. The case where k < oo can be recovered from the k — oo 
case just by relaxing the no-crossings condition. In the opposite limit (k = 0) crossing configurations 
occur freely. We will argue later how to extract conclusions for the whole range of k. 

In figure El we see loop configurations without and with intersections. The rules needed to label all 
(occupied or not) edges, and the possible vertex present in our model are shown in figure H Since our 
loops are not oriented we cannot define in a consistent way the alternated labeling of the segments of the 
loop. This is why all segments in the loop have to have the same label B. In the same way, considering the 
'dual loop' constructed from the unvisited edges, we see easily that all unvisited vertexes have the same 
label A ^ B. Now, looking at the total height difference when going around one site in our lattice, we 
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Figure 6: An example of a dense-packed non-oriented loop model configuration on a bicolorable regular 
tetravalent lattice. On left hand side there is a non-self-intersecting loop model while right hand side the 
example shows a configuration with self-crossing loops. 

can extract some consistency conditions for our height unit increments A and B. These local consistency 
conditions are depicted in figured It is easy to see then that the equations required from consistency of 
the height variables 

2A + 2B = (4) 
AA = (5) 
45 = (6) 

constrain our height variable so much that it is in fact constant all over the lattice. 

We should remember that there is no temperature or Hamiltonian in these loops models from which 
we extract the height model. The connection with our model appears when we consider the relevance of 
the bending vertex. We obtain the same equations for A and B either allowing or not the bending of 
the loops. This corresponds to studying our loops inherited from the gonihedric spin model at T = or 
T = oo, respectively ^. We should then come to the conclusion that there is no continuous transition at 
finite temperature since both extrema share the same critical behavior J^. The same reasoning with the 
self-crossing vertex being allowed or not brings us to the same conclusion for every value of the coupling 
constant k in the gonihedric Hamiltonian®. Thus we end up again in a c = theory in the whole k space. 
This means no criticality, in perfect agreement with numerical simulations [3]. 

^Note that in the loop/height model, a bending vertex has the same weight as the straight loop vertex and the bulk. This 
corresponds to infinite temperature. When we do not allow bending to appear, it is equivalent to freezing the gonihedric 
loops at T = 0. 

^Remember that k = is in fact trivially solved |in| and does not possess any thermodynamical transition. 
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Figure 7: Rules defining the height increments from the loop configuration and consistency equations 
from the different kind of vertex. In parenthesis the vertex appearing with the self-intersections. 



4 The two-dimensional gonihedric spin model coupled to gravity 

We are now going to couple the gonihedric loop model to gravity. This is interesting by itself and in 
order to understand the effect of placing the gonihedric spin model on a 'fluid' surface. Of course the 
resulting model is nothing but the one-dimensional version of the gonihedric string in a two-dimensional 
embedding curved dynamical space. ^ 

The energy of the gonihedric spin model can still be written, on a flxed quadrangulation, as 71,2 4- 4^714, 
where n2 is the number of unit 'corners' of the inter-phase between plus and minus sign spins, and is 
the number of self-crossing points of the inter-phase. 

In order to couple this model to gravity we need to put our spins in the sites of a random lattice 
built from square pieces*. We will consider, when coupling the model to gravity, the limiting case where 
the loops never cross themselves. This non intersecting limit corresponds to the k ^ cxi limit in the 




Figure 8: Example of a random lattice with a gonihedric spin model on it 

gonihedric spin action. In fig|Slwe can see an example of this kind of quadrangulations. 

To describe these configurations we shall introduce N x N Hermitian random matrices. Let us see 
how to do this. We have first of all plaquettes where no loop goes across them, then plaquettes crossed by 
one loop without bending through it (a straight loop piece across a unit square of the lattice), and finally 
plaquettes were the loop crossing them bends 'right' or 'left'. These three building blocks of the random 
lattices are presented graphically^ in figlHl with the corresponding term of the associated matrix model 

^In higher dimensional cases related studies in the continuum have been performed by 1171 . 

*We chose squares as the building blocks of our random lattice since our original model is defined in a square regular 
lattice. It happens to be important in order to define what exactly means bending the loop. 

^To simplify the visual identification with the loop model we have not drawn the lines corresponding to the bulk 
propagator, i.e. to the A matrix propagator. The loop is generated with the B matrix propagator. 
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that will generate them upon integration. As we can see we are considering the most general case where 



gABAB 
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Figure 9: Correspondence between the loop pieces and the matrix interaction that are going to generate 
them. 

all the couplings are different. In our specific case we will impose the condition g' ~ g due to the fact 
that a straight piece of loop does not contribute with any amount of energy to the action, so the coupling 
needs to be equal to the bulk coupling. In addition to these 'interaction' terms we need to include the 
kinetic term for the two matrices, i.e. the quadratic terms and that generates the propagators. So 
finally the matrix model that will represent our gonihedric loop model coupled to gravity is 

n(JV) 



ZN{g,X) = / dAdBexp 



-NTt { -(A^ + B^) + ^A'^ + ^ABAB + XA^B' 



The numerical factors preceding each interaction term are in fact symmetry factors. To our knowledge 
this model has not been solved exactly. Although somewhat similar matrix models have been solved 
elsewhere|13j. their solution cannot be applied directly to our matrix model. 

The coupling constants g, g' and A have a definite meaning in terms of the statistical model we are 
treating, g = g' are in fact representing the bulk energy of the unit square, i.e. it represents the conjugate 
variable to the area of the lattice, the cosmological constant. On the other side the combination ^ is 
equal to the statistical Boltzmann weight of each bended loop piece, that contains all the information 
from the energy of the statistical model. 

9-9'= , ^ = e-'^^-- (7) 

4.1 Partial integration 

The B matrix integration is in fact Gaussian so it can be integrated exactly. This integration was first 
presented in We first diagonalize the A matrix 

NTt ( -(A^ + B^) + ^A* + ^ABAB + XA^B"^ 



Zn{9,X) = j dAdSexp 

= / mda,dB, J [ W dReSy dlmS,, I (a, - a 

\i=l / \j>i=l / 0>t=l 



(8) 



where {oi} are the eigenvalues of A. The action S{{ai\, {Bij}) is 

^ /I \ 1 ^ 

5({a,},{i?,,}) = 5: -a2 + fan+-^B2(l + (.9 + 2A)a,2) (9) 



2 ' 4 W 2 

4=1 ^ / 1=1 

N 

+ ^ (ReB^. + ImB^.) {ga.aj + X{af + a^)) (10) 

i>i=i 

Rescaling the B matrix entries correspondingly we transform the B integral in a bunch of decoupled 
Gaussian integrals 

„ N N N 

ZN{g,X) ^ JY[da. ]^ („^._„^)2 ]-[(5^^.+g„,^„^.+^(„2^„2))-l/2g-^i:«,(ia? + fa^^ 

dA{-Det[l(g>l + 9A(g>A + XiA^ (g>l + l(g> A^)]y^ q-^^' i^^'+i^") (U) 
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where Det means that a determinant of a tensor product matrix, i.e. a iV^ x iV^ matrix. 

This is a nice expression of the two matrix model as a one matrix modeL Ahhough some especial 
cases of this model with specially fine tuned coupling constants have been solved, it appears to be very 
difficult to solve in the general parameter setting which is interesting to us. So we have to turn to other 
techniques. 



5 The renormalization group approach 



The renormalization group approach for random matrices was first introduced in |14) by Brezin and 
Zinn- Justin and later improved by Higuchi, Itoi, Nishigaki and Sakai in |19| . The improvement in jl9| 
provides a more quantitative analysis of the critical structure in parameter space. It consists basically 
in using the equations of motion of the matrix model (also called loop equations) in order to deal with 
redundant operators in the renormalized action in the large N limit. This more complete analysis is much 
more involved and we have not applied it to our model but instead used the simpler Brezin-Zinn- Justin 
approach since it is equally valid to get the qualitative flow we are looking for. 

The renormalization procedure in matrix models consists in the integration of the last row and column 
of the matrices, so that we relate the partition function zf^~^^^ to Z*-^'. One could in principle do so 
with the one-matrix reduced model 1)11(1 . but due to non linearities in Tr (A) and Tr (A^) already in the 
saddle point equation we decided to return to the original two-matrix model, where this non-linearity 
is not present. In order to make easier the 'last-row-last-column' integration we are going to write the 
matrices A and B as follows 



An+1 = 
Bn+1 = 



An 


Bn 



(12) 
(13) 



Since the action is Gaussian in the B matrix we know that the integral in the v vector and that of the 
/3 variable will be easy to perform. Only a saddle point in the a variable will be left as in the original 
approach JA^. 

To perform the renormalization group computations we need to free our parameters from the constraint 
that the original statistical model imposes on them, i.e. we should let g ^ g' . The partition function will 
satisfy 

ZN+iig,g',X) = CiNf" ZNig + 5g,g' + Sg',X + 6X) (14) 



where Sg, Sg' and SX depend on N. This equation will be well approximated by the linear differential 
equation 

\r 9 „ d ^ d ^ d 
!'dN-^^d-g-^^'W-^'dX+\ 



FiN,g,g',X)=r{N) 



(15) 



where r{N) is related to C{N) and 

FiN,g,g\X) 
is the free energy of our model. 



iV2 



log 



ZN{g,g',X) 



Zjv (0,0,0) 



(16) 



5.1 Last-column-last-row integration 

Consider 

r-(W+l) 



ZN+i{g,g',X) = 



We rotate both matrices with the same U € U(n) unitary matrix in such a way that the A matrix 
diagonahze. One unitary integral decouples and contributes with a factor Vjji^N+i) = •7r^^^+-^'/^/ Hp^^ 
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The matrix A is now 



A 



N+l 



ai 



V 



a J 



After this we rotate back the first N rows and columns so that A and B will look like Eq. (|13|l . 
It is easy to see that all the interaction terms in the action can be written as 



Ttn+i[A"] 
Trjv+i[S2] 
Ti N+i[ABAB] 



TrAr[A"] +a" 
Ttn[B'^]+ + 2v''v 
Ttn[ABAB] + 2av'^Av + a^0^ 



After substituting this into the partition function we find several terms. Those depending on v, (3, and 
a will be integrated out so that a new effective action for the N x N two-hermitian matrix model will 
arise. After performing the integrations on v and f3 we find 



ZN+i{g,g', A) 



Vr 



U{N+1) 



Vr 



C 



N 



U(N) 



N 



where 



and 



AAABe 



Cn — 



1(9,9 , A) = da 



N+l 



det[^ 



N 



(17) 



(18) 



(19) 



det[l + Aa2 + gaA + AA2](1 + (2A + g)a'^)i 

We can already see here that wc are not going to receive any direct contribution to the renormalization 
of the A or g constant in the action. There is no B matrix dependence in the saddle point integral, so 
we cannot generate the operators Tr (A'^B^) or Tr (ABAB) in our renormalization scheme. The only 
flow of this two parameters comes from the renormalization of the "kinetic" Tr {A"^) operator after field 
rescaling. The a integral cannot be computed explicitly. The saddle point method has been used to solve 
it in the large N limit. The two determinants in the integrand contribute to the saddle point while the 
square root in the denominator does not. 



5.2 Saddle point equation 

The saddle point equation reads 

N \as-Aj V ^ TV \l + Xa^^+gasA + XA^J ^ ' 

To solve this saddle point equation we will assume that a^, the solution of the saddle point equation, 
has an expansion in powers of aj = Tr [A^). Once this saddle point is computed we must introduce 
as into the integral, then we will again expand everything to find the contribution of this integral to the 
operator 0,2 

I{9,9'A) = dei[A-asf r exp -(A^ + 1) f ^ + ^a^) (21) 

dct[l + Xal+gasA + \A^]{l + {2\ + g)al)^ ' 4 7 

= exp-(A^ + l)(Co + Ciai + Ci,ia? + C2a2 + ---) (22) 

To keep track of the terms we need in the expressions we will perform the change of variables A eA. 
At the end of the calculation we shall put e = 1 . 

nr.^ f 1 \ , I z\ , ^ nr. f 2\as+geA 



— Tr \ ={as+ g'al) + — Tr ^ — ^ . (23) 
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We are going to expand everything up to second order in e. This will be enough to extract the C2 
coefficient. We expand 

as^af^ +ea^P +^af'> +0{e^), (24) 

where 

5(1) = aWai, 5(2) = af^a^ + a^^-'^^a\. (25) 
At order ©(e*^) we get an equation for ah^ 

<^f^ + - \ ro.. = (26) 

a 



(0) ^ ;,^(0)3 



that is a sixth order equation reducible to a third order one. Of all solutions we will chose the one that 
matches the known solution for g ^ and A ^ 0, that is ai°^ — > \/2 and which is real. At next order 
we get a linear equations for ai^\ and at order 0{e^) we get two linear equations. We have tested the 
stability of the saddle point solution we have found in the vicinity of the point g = X = using the above 
expansion. We have checked that the solution is stable for small values of these coupling constants in the 
whole g' range. 

At this point we can express the saddle point evaluation of the integral in terms of these constants 
and find the renormalization of the operator Tr (A^) = Na2- If we call Z2 the coupling relative to that 
operator (when we begin Z2 = 1) the renormalized coupling Z^*^"^ after one renormalization group step 
from the original model will be 

^ren ^ ^^±1 f ^ _^ _}_C^ (27) 

where 



2 + %al - (ff^ - A X^)al 
2(ao + Xal) 

and ao is the solution we took from the sixth order equation Eq. H26|) . 



'-2 - — , x„3A2 (.28) 



5.3 Renormalization-group flow 

We must scale the A matrix in order to keep Zj*^" = 1 at the end of the complete renormalization group 
process. This induces a flow on the couplings A and g. 

= A (^1 - + 0(Ar-2)^ (29) 

ff'-™ = g (^1 - 1^2 + 0(7V-2)) (30) 

The behavior under the renormalization-group depends on the sign of Ci ■ 

If we pay attention to the statistical meaning of the couplings in the original model, it is clear that 
we can keep g ^ 1 and both g and Xj g are positive, i.e. X,g > is the physical quadrant. As in the 
usual renormalization group approach ^^glg, ^ < near g = g* implies an infrared fixed point a,t g ~ g*. 
This is the case in the physical quadrant. 

In other words, the two constants A and g renormalize towards its fixed value X = g = 0. Thus only 
the g' coupling constant remains. The constant g' accompanies is the operator Tr {A^) which generates 
only pure gravity critical behavior. This means Cmatter = 0. 

We can also think of this flow in the following way. Regardless of g, if A is flowing towards zero, we 
can forget about the corresponding operator and concentrate on the rest of the action. The corresponding 
model has been solved in '20' . In our special choice of parameters g' — g this brings us to pure gravity 
again. 

The irrelevance of the operator Tr {A^B^) could have also been already predicted with the heuristic 
tools presented in the first half part of this paper. Using the height model language from we can see 

^"We could also choose — V2. 
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Figure 10: and ^ for positives values of the coupling constants around the origin (physical quadrant). 
Two different values of g' have been plotted to show the positiveness of this quantity for all g'. (a) g' ^ 0. 



that for ordinary gravity'^-'^ the height model we work with contains exactly the same constraints whether 
we include the Tr (A^B^) operator or not. These arguments show that this operator is not modifying 
the essential nature of the scalar model living in the continuum. Thus the central charge should be the 
same. 

Although a complete solution of this matrix model has not been achieved, we have been able to explore 
the renormalization group flow. We find again the conclusion that the gonihedric loop model coupled to 
gravity corresponds to c = central charge for the matter part^, i.e. pure gravity. This again agrees 
with the lack of a continuous thermodynamical transition observed for the two-dimensional gonihedric 
spin model. We therefore find a rather nice agreement among all approaches. 

6 Conclusions 

We have presented two, quite different, versions of the gonihedric loop model. 

First, through a direct duality we have transformed our gonihedric spin model into a loop model. We 
have studied the gonihedric loop model through the height models using the results given in ^1 . 
This height models has been proved to be give very satisfactory predictions for the central charge and 
critical exponents for combinatorial loop models. 

We have applied height model technology to the present problem and we have extracted the central 
charge in some limiting cases of the gonihedric loop model. These limiting cases correspond to zero and 
infinity temperature, and zero and infinity k (self-avoidance parameter). Modulo some very plausible 
hypothesis concerning the interpretation of the central charge and renormalization-group flows, we are 
then able to conjecture the value for the central charge found for the limiting models to the whole 
temperature and k space. The value of the central charge happens to be the same for the whole space 
of parameters (c = 0), so we find no continuous transition for the gonihedric model (in agreement with 
existing numerical simulations 9 ). 

Secondly, we have coupled our loop model to gravity in the k — > 00 limit. We have worked out a 
partial integration since in this limit the loop matrix is Gaussian. Unfortunately we could no go further 
in the exact integration of this matrix model due to its strong non- locality (in eigenvalue space). 

A renormalization group approach has proved useful to bypass this problem. We have computed the 

^^respect Eulerian gravity. See 1121 . 



(b) g' = 1. 
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infinitesimal renormalization flow for the two couplings g, and A. We found that they flow always towards 
zero, so only a quartic one-matrix model remains. We can also consider only A as irrelevant coupling 
constant. The limiting matrix model has been already solved and it turns out to correspond to pure 
gravity. Then using the KPZT^ relations we rediscover the fact that c — Q for the two-dimensional 
gonihedric spin model. 
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